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ABSTRACT 

Most models of volatile delivery to accreting terrestrial planets assume that the carriers for water 
are similar in water content to the carbonaceous chondrites in our Solar System. Here we suggest 
that the water content of primitive bodies in many planetary systems may actually be much higher, 
as carbonaceous chondrites have lost some of their original water due to heating from short-lived 
radioisotopes that drove parent body alteration. Using N-body simulations, we explore how plan¬ 
etary accretion would be different if bodies beyond the water line contained a water mass fraction 
consistent with chemical equilibrium calculations, and more similar to comets, as opposed to the more 
traditional water-depleted values. We apply this model to consider planet formation around stars of 
different masses and identify trends in the properties of Habitable Zone planets and planetary system 
architecture which could be tested by ongoing exoplanet census data collection. Comparison of such 
data with the model predicted trends will serve to evaluate how well the N-body simulations and the 
initial conditions used in studies of planetary accretion can be used to understand this stage of planet 
formation. 


1. INTRODUCTION 

Low mass stars are attractive targets for the search 
for habitable worlds due to their great abunda nce in the 
galaxy (nearly 75% of stars are M dwarfs, iScalo et all 
1200711 . and the fact that planetary detection signals from 
radial velocity, transits, direct imaging and microlens- 
ing would be greater for these stars than those with 
higher mass and greater size. As such, recent efforts 
have been dedicated to predicting the characteristics of 
planets that may form around such stars in order to eval¬ 
uate how common habitable planets are in the galaxy 
(e.g. Ray mond et al.ll2007t lLissaue dl2007t [Qgihara fc Idal 
2009t iMontgomerv fc Laughlinll2009lh 

One of the primary objectives in these investiga¬ 
tions has been to evaluate the likelihood that plan¬ 
ets around these stars may have accreted water dur¬ 
ing their formation and th us b e potential abo des for 
life. [Raymond et al.1 (|2007f l and iLissaued (120071) exam¬ 
ined this possibility in detail, considering whether the 
same processes which we believe were responsible for the 
delivery of water to the young Earth could have also oper¬ 
ated around these lower mass stars. That is, Earth is ex¬ 
pected to have accreted largely from native, dry planetes- 
imals and embryos as temperatures were too high at 1 AU 
in the solar nebula for water to have condensed as a solid. 
Earth’s oceans would thus have come from water-bearing 
solids which formed outside the water line-the distance 
from the young Sun where water ice was able to condense 
and be incorporated into solid planetesimals. The accre¬ 
tion of such materials from more distant regions is a natu¬ 
ral consequence of the dynamics of planetary accretion in 
our Solar System (|Morbidelli et all 120001 lO’Brien et al.1 


120061 : iRavmond et al.l [20 09). and is consistent with the 
similar D/H ratio observed in Earth’s oceans and clion- 
dritic meteorites whose pa rent bodies now populate this 
regio n of the Solar System (|Alexander et al.ir2012t iMartvl 
120121 ). _ 

IRavmond et al.l (1200711 first explored how much water¬ 
bearing material was accreted by habitable zone (HZ) 
planets that formed around stars of different masses. 
Scaling the location of the wat er line with th e main 
sequence luminosity of the star, IRavmond et al.l (120071) 
found that while the frequent gravitational encounters 
between bodies beyond the water line around solar-mass 
stars would lead to some of these materials being scat¬ 
tered inward to be accreted by the inner planets, such en¬ 
counters were less frequent and less intense during planet 
formation around low-mass stars. As a result, little to no 
water-bearing bodies would be scattered inwards, leaving 
the planets in the habitable zones around these stars to 
accrete predominately native, dry materials, thus leaving 
t hem vo latile poo r. 

iLissauer (2007 ) pointed out an additional difficulty for 
the accretion of volatiles by HZ planets around low-mass 
stars. Specifically, the slow contraction of low mass stars 
means that they would be more luminous during the 
epoch of planet formation than during their main se¬ 
quence lifetime. This increased luminosity is important 
as irradiative heating by the central star plays an impor¬ 
tant role in setting the location of the water line in a 
protoplanetary disk, particularly when accretional heat¬ 
ing th rough the midpl a ne of the disk is negligible. Thus, 
while IRavmond et al.l (j2007l ) scaled the location of the 
water line with the main sequence luminosity of the star, 
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iLissauerl (12007ft argued that the water line would be lo¬ 
cated even further from low-mass stars than this scaling 
would suggest. Thus water-bearing solids would have to 
be scattered inward from greater distances in order to 
be accreted by HZ pla nets around ve r y low mass stars. 
Given the findings of iRavmond et al] (2007 ). such scat¬ 
tering would still have been minimal given the low surface 
density of solids expected in the disk, again leaving HZ 
planets water-poor. 

An important assumption in these models is that the 
water-bearing planetesimals and embryos beyond the wa¬ 
ter line in all of these systems were similar in com¬ 
position to carbonaceous chondrites, the water-bearing 
primitive bodies which are thought to be the source 
of Earth’s water dMorbidelli et at JgOOO yO’Brien et al.1 
1200111 : IRavmond et al.l 120091: lAlexander et, al.l 1 20 bi b In 
our Solar System, these bodies formed through the ac¬ 
cretion of a mix of rocky minerals and water ice be¬ 
yond the water line. Once they accreted, the decay 
of short-lived radionuclides (SLRs), most notably 26 Al, 
raised the internal temperatures of the bodies, melting 
the water ice and allowing it to react with the rock 
to form hy drated minera ls dGrimm fc McSweenl fl 989: 
iCohen fc Coked [20001 : iCastillo-Rogez fc McGordl 2010b 
This process was also capable of mobilizing or driving 
off water, meaning that the final water content seen in 
rocky meteorites is less than the original accreted mix. In 
fact, under equilibrium conditions in the solar nebula, the 
water-to-rock ma ss ratio beyon d the water line should 
be nearly 1-to-l dLoddersI 1200.3 ;: however carbonaceous 
chondrites contain only up to 5-10% water by mass, sug¬ 
gesting possible loss of this water after accretion. Thus 
in other planetary systems where 26 Al may have been 
less abundant or absent, bodies beyond the water line 
could have been much more water-rich and the total sur¬ 
face density of solids greater than previously considered 
in planetary accretion studies. 

Here we explore how a greater water content of distant 
bodies would affect planetary accretion around stars of 
different masses . Sp ecifical ly, we bu i ld on the studies of 
IR avmond et al.l (| 200 7ft and ILissauerl (|2007ft and focus on 
how the absence of live SLRs during planet formation, 
which could be the norm for many planetary systems, 
affects the accretion of planets in such systems and the fi¬ 
nal water content of potentially habitable planets around 
low mass stars. In the next section we discuss in greater 
detail the water content of planetesimals in protoplane¬ 
tary disks and the location where water is expected to 
be incorporated into solid bodies. We follow that by pre¬ 
senting a series of simulations which compare how plan¬ 
etary accretion occurs under the standard assumption of 
using carbonaceous chondrites as analogs for the bodies 
originating beyond the water line and with how it occurs 
when those bodies are composed of equal parts water and 
rock. We end by discussing our results and highlighting 
areas in need of future work. 

2. WATER DISTRIBUTION IN PLANETESIMALS 

Water condenses in a protoplanetary disk where its 
partial pressure exceeds the saturation pressure, P^o > 
Pgo- While the exact temperature where this occurs 
depends on the details of the physical structure of the 
disk and relative abundances of elements, the expected 
temperature ranges from 140-180 K for typical disk con¬ 


ditions. Here, we will adopt a value of T con d =140 K to 
guide our calculations, recognizing that the results of our 
study are not very sensitive to the exact choice of conden¬ 
sation temperature. Further, choosing a value at the low 
end of the expected range means that the results pre¬ 
sented here are conservative, with higher temperatures 
meaning the water line would be located closer to the 
star, making it easier for water-containing bodies to be 
delivered to HZ planets. 

The temperature structure of our protoplanetary disk, 
the solar nebula, throughout its evolution is uncertain, 
and it likely evolved with time as the physical properties 
of the Sun and nebula itself evolved. Nonetheless, the 
water line of the solar nebula has typically been thought 
to have left its mark at ~2.5 AU, as C-type asteroids, 
which are thought to be the parent bodies of carbona¬ 
ceous, water-bearing meteorites, are the dominant type 
of meteorite found outside this distance, while dry, S- 
type aste roids, are the main bodies found interior to this 
distance dGradie fe Tedescol 1982 : rMorbidelli et al.ll2000l : 


IRavmond et al.l 120051120071 2006). However, it is worth 
noting that recent observations indicate that this zon¬ 
ing is only true for the largest objects, and that asteroid 
compositions exhibit greater variations when smaller ob¬ 
jects are considered (|DeMeo fe Carrvl 12014ft . This has 
been interpreted as evidence that the positions of the 
asteroids today are not the same as their formation loca¬ 
tion, and instead the asteroids were scattered into their 
current positions as the planets were assembled, some¬ 


thing that h as been suggested bv a variety of dynam ical 
models fe.g. iBottke et al.ll2006at iWalsh et aI1l2012ft . If 
so, then the distribution of asteroid compositions today 
cannot be used to provide strong constraints on the ther¬ 
mal structure of the solar nebula. 

In practice, the temperature within a protoplanetary 
disk, and thus the location of the water line, is deter¬ 
mined by balancing the heat inputs from irradiation from 
the central star, the ambient radiation field, and viscous 
dissipation within the disk with the energy lost by ra¬ 
diation at the disk surface. A key issue in determining 
the thermal evolution of the disk is identifying the mag¬ 
nitude and location where viscous dissipation occurs. If 
dissipation occurs throughout the entire thickness of the 
disk, temperatures everywhere will be increased because 
of the added energy. However, it is unclear whether dis¬ 
sipation would occur throughout, or if it would be lim¬ 
ited to particular regions of the disk. For example, if 
the MRI is responsible for driving disk accretion, vis¬ 
cous heating may be limited to regions of high ionization 
fractions like the extreme inner edge of th e disk, t he disk 
surface layers, or the outer disk (e.g. lGammielll996ft . In 
such a case, irradiation from the central star is likely the 
dominant source of heating in the regions where planets 
formed, and thus would set the location of the water line 
(|Lecar et alJl2006l : iMartin fe Livioll2Q12ft . 

When irradiative heating is dominant, the midplane 
temperature at some distance, r , fr om t he cen tral star 
can be calculated as (jGhiang fc Goldreichl 1 997ft : 



(1) 


where T* is the effective temperature of the star, i?* is 
the radius of the star, and a is the grazing angle of the ra- 
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diation on the disk surface. The grazing angle will likely 
depend on location in the disk, as it is determined by 
the height of the dust layers in the disk which will evolve 
with time as dust grows and settles to the midplane. De¬ 
spite these details, a is expected to have a value of ^0.05 
in the regions of the protop la netar y disk where the water 
line is expected (iChiang fe Goldreichl 119971) . with slight 
differences being unimportant as temperature is propor¬ 
tional to a 4 -an angle that is 50% larger would lead to 
a temperature change of just ~T5%. Again, if viscous 
dissipation within the disk were taking place throughout 
the disk, then temperatures would be greater than this, 
and thus the temperatures given by this equation should 
be considered a minimum value. 

Figure 1 illustrates how the location of the water line 
evolves in a disk as a result of the pre-main sequence 
contraction of stars o f dif feren t masses. Here the stellar 
tracks from ISiess et al.l (120001) . which calculate the evo¬ 
lution of the radius and temperature of stars of different 
masses over time, were used along with Equation 1 to 
determine the region in the disk where T=140 K. Note 
that other stellar models would likely give slightly dif¬ 
ferent numbers, though the trend of cooling with time 
are robust. The change in water line location is most 
rapid early in stellar evolution as the star contracts at 
the greatest rate. 

As planetary accretion simulations determine the 
provenance of the different building blocks for a planet 
in order to determine its final composition, the location 
of the water line marks a transition from dry planetesi- 
mals and embryos to those that are water-bearing. Pre¬ 
vious models for the accretion of the terrestrial planets 
have largely used carbonaceous chondrites as analogs for 
the bodies originating just beyond the water line and 
thus have assigned a water conten t of ~5% to those bod¬ 
ies coming from that region (e.g. IMorbidelli et al.l f2000t 
IQ’Brien et aD B IRavmond et all 120051 120071 120090 . 
Indeed, carbonaceous chondrites contain 5-10% water by 
mass, though this water is not stored as pure water or 
ice, but instead within hydrated silicates. These miner¬ 
als formed on the meteorite parent bodies, as the de¬ 
cay of 26 A1 heated the parent body, melting any wa¬ 
ter ice which was present (melting temperature ~ 273 
K, thus >100 K higher than the nebular condensation 
temperature). The liquid water then reacted with the 
anhydrous rock, forming minerals such as serpe ntine, 
brucite, and cronstedt i te (e.g.lGrimm fe McSweenl 1989; 
Cohen fc Coked 120001 ICastillo-Rogez fc McCordl 2010 : 
Hcw ins et al.ll20l4 lZolotovll2014D . 

Thus it is important to remember that the ~5% mass 
fraction of carbonaceous chondrites we see today is the 
product of the geophysical and geochemical evolution of 
the original parent body which formed from a mix of 
rock and ice over 4.5 billion years ago. It is likely that 
the original parent bodies contained more water than 
seen in the carbonaceous chondrites today. The same 
radiogenic heat which melted the water ice and allowed 
hydrated minerals to form also could have mobilized the 
water, allowing it to migrate to the surface of the aster¬ 
oid where it was lost to space or frozen out in the surface 
layers where it would later sublimate or be eroded by 


impacts ( Grimm & McSween 1989k I Cohen & Coked 

2000; 

Young et al. 20031: Castillo-Rogez & McCord 2010 

). In 


fact, cronsted tite , which is found in many hydrated me¬ 
teorites (e.g. iHewins et al.l I2014D is thought to require 
very high water-to-rock ratios (wa ter-to-ro ck mass ra¬ 
tio >1-10, or even higher) to form (lZolotovll2014 ). Fur¬ 
ther, based o n the mineralogy of CO carbonaceous chon¬ 
drites, [HowarT£r^I3 ( 2014 ) has estimated that these me¬ 
teorites formed in bodies that had a water mass fraction 
of ~50%, with water being lost through oxidation reac¬ 
tions and degassing that occurred during the early ther¬ 
mal evolution of the parent asteroid. Such mass fractions 
are not unheard of in asteroidal bodies, as the largest 
body in the asteroid belt, Ceres , is believed to contain 
at least 20-30% water by mass ([McCord fc Sotinl 120051 : 
ICastillo-Rogez fc McCordll2010D . Much of this water is 
expected to be found as the outer mantle of the asteroid, 
which for med as a re sult o f different iation of the aster¬ 
oid (ICastillo-Rogez fc McCordl l201(j ). The outer layers 
would have been subjected to loss of water by sublima¬ 
tion and impact erosion, suggesting that the water seen 
in Ceres today is only a lower limit on the amount which 
was present when the asteroid first formed. 

These high water mass fractions are consistent with 
chemical equilibrium models for the composition of solids 
in the solar nebula which predict that solids bey ond 
the w ater line should be ~50% water by mass dLoddersI 
120031) . Indeed, such ratios are seen in other bod¬ 
ies in the outer solar system, with comets and icy 
satellites having water mass fract i ons o f ~40-50% (e.g. 
Delsemm el 11992 : IMorbidelli et all 120001 : iSchnbert et, all 


20101 : lA’Hearnl I2011L and references within). How¬ 


ever, even higher values are possible as this estimate 
ignores any dynamic effects which could occur within 
a protoplanetary disk. Diffusion of water vapor in 
the protoplanetary disk and inward drift of icy solids 


the water line (Stevenson 

& Li 

nine 19881 ICvr et al. 

19981 SteDinski & Vala 

seas 

1997T 

Hueso & Guillot 20031 

Cuzzi & Zahnle 20041; Ciesla & Cuzzi 

2001). The en- 


hancement of water in this region could lead to even 
greater water mass fractions than expected under what 
is predicted under chemical equilibrium. 

Thus, here we explore how water-rich planetesimals 
from immediately beyond the water line would impact 
the accretion and chemical properties of the terrestrial 
planets that form around a star. Water loss from plan¬ 
etesimals in our solar system was likely due to radiogenic 
heating, notably the decay of 26 Al, whose abundance 
was likely variab l e acro ss different planetary systems. 
iVasileiadis et al.l (I2013D demonstrated that the 26 Al 
abundance in giant molecular cloud cores will increase 
in time as early formed stars introduce freshly synthe¬ 
sized short-lived radionuclides to the star-forming region. 
Thus early-formed planetary systems would be lacking 
the radionuclides to drive thermal evolution on planetes¬ 
imals, while those that formed later could be similar to 
our own Solar System, allowing s imilar thermal evolution 
of sm all bodies to take place (Uura et al.l 120131 : lYoungl 
I2014D . However, the efficiency with which radionuclides 
would be injected and mixed into a forming planetary 
system is uncertain , and could be rare or require special 
circumstances (e.g. IBoss fc K ciscr 2 0Tot lOuellette et al.1 
120101 : iWilliams fc Gaidod 120071 : lAdams et al.l I2014D . As 
a result, it is necessary to consider how the absence of 
such radiogenic heating could preserve water-rich plane- 
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tary building blocks beyond the snow line, making them 
more comet-like in their water-to-rock ratios than car¬ 
bonaceous chondrite-like, and the consequences for plan¬ 
etary accretion. 

3. PLANETARY ACCRETION SIMULATIONS 


In order to explore the effects that water-rich planetesi- 
mals would have on planetary accretion, we perfor med 40 
simul ations using the MERCURY N-body code dChambersI 
EMI). As the primary goal here is to understand how 
the presence of water-rich planetesimals would affect the 
volatile inventory of planets around different stars, we 
follow a similar approach as iRavmond et al.l (120071 ) to 
facilitate comparisons between our findings and previous 
results. That is, we consider stars ranging in mass from 
Af»=0.2-1.0 Mq and assume that the surface density of 
rocky planetary building blocks is described by: 


£ (r) = £i 



where £i is the surface density at 1 A U, and is the 
mass of the star of interest. Following IRavmond et al.1 
(j2?m we adopt a value of £i =6 g cm 2 for all cases con¬ 
sidered. Note that we are assuming here that the metal- 
licity and relative elemental abundances of the plane¬ 
tary systems are identical to those of our Solar System- 
t he impact of varia tions in elemental abundances (e.g. 

I Johnson et al.l 12013 1 on the properties of planets is left 
to future studies. For each system, we assume the solids 
in the disk are distributed from an inner radius, Ri n , 
to an outer radius R ou t- While the choice of such pa¬ 
rameters are largely arbitrary in these models, we scale 
Rout to correspond to the location where £ {R ou t)= 0-75 
g cm -2 . R in was chosen to be roughly half the distance 
inward of the habitable zone. T his differs from the ap¬ 
proach of IRavmond et al.1 (|2007f ) who scaled the inner 
and outer boundaries of the solids disk with the main 
sequence luminosity of the star, though the inner bound¬ 
aries used here and in that study are nearly identical, 
while the outer edges of the disks are generally at greater 

distances in our simulatio ns. _ 

Where IRavmond et al.l (l2007t) considered all of the 
solid mass in the disk to be contained within plane¬ 
tary embryos, here we distribute 50% of this mass into 
embryos and 50% into planetesimals to represen t the 
initial stages of oligarchic growth (iKokubo fc Rial IT9981 
lO’Brien et al.l 120061 : IRavmond et al.l 120091) . As" in pre¬ 
vious studies, embryos interact gravitationally and col- 
lisionally with both other embryos and planetesimals, 
while planetesimals only interact with the embryos. We 
assume that the initial mass of the embryo is given by 

M em b=Aj Af$; that is, for a solar mass star, the 

embryos have a mass roughly half of that of Mars, with 
this mass decreasing proportionally with mass of the star. 
In all cases, the mass of planetesimals in our simulation 
is equal to ^ that of the embryos. Previous models have 
varied in how they have treated the initial masses of the 
embryos: some use variable embryo ma sses spaced by a 
multiple of their mutual Hill radii (e.g. IRavmond et al.l 
I2007ll2009h while others have used embryos of constant 
mass spaced in a way to mirror the surface density of 
the s olids of the disk (e.g. lO’Brien et al.l [20061 : ILissaueil 
2007). While details about individual runs are known 


to vary with the exact initial conditions of an N-body 
run, the trends in regards to the final outcomes (planet 
number, masses, volatile content) that emerge when ex¬ 
amining these different studies are independent of the 
choice of how embryo properties are initialized. Thus we 
expect the results described here to be independent of 
the choice of embryo mass distribution used: constant 
embryo mass or variable embryo mass. 

We consider two scenarios for each star. First, we fol¬ 
low previous studies in assuming that the building blocks 
in our models are made predominately out of rock, and 
that beyond the water line the bodies are composed of 
5% water. We label these simulations as CC (for ‘car¬ 
bonaceous chondrite’). We also consider situations where 
the water content of the materials originating beyond the 
water line is 50%, and label these simulations as IC (for 
‘icy’). In these later cases, the surface density of materi¬ 
als beyond the water line, (r > Twl) is doubled to reflect 
the increase in the mass of available solids there. 

In all cases, the water line is assumed to occur at 
T= 140 K, with its position is calculated using Equation 
(1) and adopting the ste llar parameters for 1 Myr as 
given by ISiess et ahl (j2000f ). The location of the water 
line, along with other model parameters, in each case is 
given in Table 1. Note that the relative locations of the 
water line in the different mass stars remain nearly con¬ 
stant after 1 Myr, except for the lowest mass star. The 
water line sets the boundary between the ’dry’ and ’wet’ 
planetesimals and embryos in each simulation, as these 
bodies are expected to reflect the composition of their lo¬ 
cal environment when they form within the gaseous disk 
given that the planetesimals grow from the dust present 
in the disk and embryos grow from the planetesimals. 
Note that the water line is expected to be much closer to 
the central star i n the higher mass cases than was used in 
IRavmond et al.l (|2007l h This is because the temperature 
profile of the disk is determined by the reprocessing of ra¬ 
diation at the disk surface downward. For the lower mass 
stars, h owever, the wate r line i s located much further out 
than in IRavmond et ahl ( 2007 ) as the luminosities of the 
stars during their early evolution are much greater than 
their main sequence luminosities, even offsetting the indi¬ 
rect radiation geometry caused by the disk, as discussed 
by ILissaueil (l2007f) . 

The choice of setting the location of the snow line at 
1 Myr is somewhat arbitrary. As the stars continue to 
contract with time, the location of the water line will mi¬ 
grate inwards under the assumptions used here, as illus¬ 
trated in Figure 1. In our own Solar S ystem , there is evi- 
dence that planetesim al dKruiier et al.l[2014ll and embryo 
(jDauphas fe Pourmandl 120111) formation took place <1 
Myr into its evolution, while some meteorites prov ide evi¬ 
dence of parent body form ation 2-5 Myr later (iKrot, et al.1 
12003 IKleine et aU feonAV Formation of Jupiter via core 
accretion likely requires relatively rapid planetesimal for¬ 
mation given that envelope accretion takes millions of 
years, possibly lon ger than the median lifetime of pro¬ 
toplanetary disks (Hubick vi et al.l [2005). Thus plan¬ 
etesimal and embryo formation likely spanned millions 
of years within the solar nebula, meaning the bound¬ 
ary between dry and wet planetesimals was likely not 
static, but evolved with time, and this extended period 
may have impacted th e chemistry of the forming planets 
(jMoriartv et al.lfeOlJfl . Here we follow previous studies 
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by using a singular location for the water line to represent 
an average location. As discussed here and in previous 
studies, these results should be used to be a general guide 
on how planetary accretion would occur under the con¬ 
ditions described here as every system likely varied in 
terms of the stages of evolution that preceded the stage 
studied here. 

4. RESULTS 

Figure 2 shows the dynamical evolution of the em¬ 
bryos and planetesimals for one of our O.4M 0 CC cases, 
while Figure 3 shows a simulation for one of the IC cases 
around a star of the same mass. Accretion in each of 
these simulations proceeds as has been documented in 
previous studies: the orbits of planetesimals and em¬ 
bryos are excited by gravitational interactions with one 
another leading objects to be put on more eccentric, and 
thus crossing, orbits. The crossing of the orbits leads to 
the formation of larger bodies as a result of the collisions 
that occur. In general, planetesimals achieve more ex¬ 
cited orbits as seen by the higher eccentricities of these 
smaller bodies in Figures 2 and 3. This dynamical fric¬ 
tion has been seen in previous studies (e.g. lO’Brien et all 
120061 ) and means that these small bodies are much more 
likely to be accreted by planets which formed at differ¬ 
ent semi-major axes than the planetesimals. Thus these 
planetesimals are likely to carry materials across greater 
radial distances than the embryos, and deliver materials 
formed at more distal locations than the planets. 

This last point is important in considering the differ¬ 
ences in the planetary systems which result in the final 
panels of Figures 2 and 3. In the CC case, the inner most 
planets have negligible water content as they accreted 
largely from native materials, that is materials that orig¬ 
inated in the same region where those planets reside at 
the end of the simulation. Even with the inclusion of 
planetesim als here, this result is in agreement with the 
findings of lRavmond et all (I2007T) . who also found that 
planets which formed close to the star around O.4M 0 
mass stars would largely be composed of dry materials 
from the inner regions of the solid disk. This contrasts 
with the results of the IC case, where the inner most 
planets have water mass fractions of 1.2-3.3%. As dis¬ 
cussed further below, the ability for inner planets to ac¬ 
crete water-bearing planetesimals is due to the greater 
mass of material present beyond the water line in the IC 
cases compared to the CC cases. In the case of the two 
inner most planets, their water content was reached by 
accreting one single planetesimal each from beyond the 
water line in addition to the native materials it accreted 
during their growth. The third planet in the system ac¬ 
tually accreted 6 planetesimals from beyond the snow 
line, leading to a greater water abundance. Unlike the 
first two planets, the third planet was also formed by the 
merger of two embryos, those being the third and fourth 
embryos when ranked in terms of initial semi-major axis 
at the beginning of the simulation. These embryos were 
both dry, having originated inside the water line, mean¬ 
ing that while the third planet accreted 6x the mass 
of water-rich planetesimals, the fraction of water in the 
planet was only 3x that of the inner most planets. 

These trends continue when we consider the full suite 
of simulations around stars of different masses. Figure 
4 shows the final planetary systems that were formed 


in each of the cases considered here, with stars of 0.2, 
0.4, 0.6, 0.8, and 1.0 M 0 . While some planetesimals 
remain unaccreted at the end of each simulation, and 
thus could potentially be accreted at later times (much 
like the planets in our Solar System today continue to ac¬ 
crete extraterrestrial materials), these planetesimals have 
largely diffused to larger semi-major axes and thus are 
less likely to have a significant impact on the planets 
which formed close to the star, in the habitable zone re¬ 
gion. In each case shown, the four CC simulations are 
shown in the bottom of the panels, while the IC cases are 
shown above them. The Solar System is shown at the top 
of each plot for reference, with the water mass fractions 
for the planets set at 10~ 5 for Mercury, 3.7xl0~ 5 for 
Venus (|Donahue et al.l[T997f) . ~0.001 for Earth (a lower 
limit a s more water may be contained in the Earth’s 
mantle iM orbide lli et all fcoQOf) . and 3.5xl0~ 4 for Mars 
(jLunine et all l2003lh Again, we stress that the Solar 
System is provided for reference-our goal here is not to 
reproduce the properties of our own Solar System, but 
rather to explore the diversity of planetary systems that 
may develop and how this is impacted by the enhance¬ 
ment of water ice beyond the water line in different plan¬ 
etary systems. 

Among the key differences that are seen in these simu¬ 
lations are that the water content of the innermost plan¬ 
ets around each star are greater in the IC cases when 
compared to the traditional CC cases. These wetter plan¬ 
ets are due to the increase in mass beyond the water line 
in the CC cases. Firstly, the fact that the water con¬ 
tent of each planetesimal in the IC case is 10 x greater 
than the traditional water content of planetesimals con¬ 
sidered in the CC case means the accretion of just one 
single body from beyond the water line would deliver sig¬ 
nificantly more water than considered in previous work. 
Further, the extra mass beyond the water line (in the IC 
case there is 2 x more mass than in the CC case) results 
in more frequent and stronger gravitational interactions 
among solids in this region, leading to more significant 
scattering of planetesimals. This allows such bodies to 
undergo much more radial mixing than in the CC cases. 
We note that as in previous studies, we are assuming per¬ 
fect mergers in these simulations and not allowing for the 
potential water loss that could occur during accretion. 

It is worth noting that while the CC cases considered 
here result in drier HZ planets than the IC cases, these 
planets t end to contain more vo latiles than found in the 
study bv iR.avmond et all ( 2007 1. In particular, here the 
HZ planets that form around the O.6M 0 stars contain 
volatiles in the CC cases, while HZ planets which f ormed 
aroun d the O.8M 0 stars were fairly dry in lRavmond et all 
lj2007f ). This difference is due to the initial conditions 
of these simulations including 50% of the solid mass in 
planetesimals, which were neglected in the previous stud¬ 
ies. This equal-mass distribution of embryos and plan¬ 
etesimals are consistent with more recent simulations of 
planetary accretion fe.g. lRavmond et aLli2009l : fChambersl 
l20l3 ). As the water-rich bodies in the previous studies 
tended to be much more massive than considered here, 
they did not experience significant radial migration due 
to gravitational scattering events. Here, however, the low 
mass planetesimals are easily excited to high eccentric¬ 
ity orbits, as shown in Figures 2 and 3, allowing water¬ 
bearing bodies to be accreted by forming planets over a 
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wide range of semi-major axes. Thus planetesimals, or 
smaller, low mass solids from beyond the water-line are 
potentially the critical carriers of water to HZ planets 
around M stars. 

A second important trend that can be seen in these 
simulations is that the number of planets that form 
in these simulations increases around the smaller mass 
stars. This t r end w as also seen in the results of 
1 Raymond et al.l ()2007t ). and arises due to the increased 
available mass at the beginning of the simulations and 
the greater level of gravitational interactions around the 
higher mass stars. That is, as described above, plane¬ 
tary accretion arises due to the orbital excitation of their 
building blocks which leads to collisions. This excitation 
occurs much more readily in the higher mass disks due 
to the greater masses of bodies involved. This also leads 
to the greater number of planets around low mass stars, 
as the lower masses in these disks do not lead to signifi¬ 
cant orbital excitation of the embryos and planetesimals. 
This trend is also seen in the Kepler data as there have 
been greater numbers of planets found a roun d individ- 
ual low-mass stars than h igh mass stars (jHoward et al.1 
l2012t IMulders et aI1l2014H . In addition, the orbital spac¬ 
ing between the planets in the low mass star cases is less 
due to the lower levels of gravitational interactions. This 
trend is also seen in comparing our CC cases to the IC 
cases-there are a greater number of planets formed in the 
CC cases due to the lower masses of the disks in these 
simulations, while the increased mass beyond the water 
line in the IC cases allows for more orbital excitation of 
the planetary building blocks. 

5. DISCUSSION AND CONCLUSIONS 

Here we have shown that the increase in solid surface 
density in a forming protoplanetary system due to the 
condensation of water ice can have an important effect 
on the properties of the planets which form around a 
young star. The increase in solid surface density leads 
to more massive planets and planets with greater water 
contents than would be found in systems where water¬ 
bearing planetesimals are treated as carbonaceous chon- 
drite analogs. These effects are not limited to the regions 
immediately outside the water line-gravitational scatter¬ 
ing means that the supply of planetesimals to regions 
inside the water line, including the eventual habitable 
zone, will be increased. This allows for more massive 
and volatile-rich planets to form closer to the star than 
previously realized. 

The findings here are generally in agreement with pre¬ 
vious studies of the volatile con tent of hab itable zone 
planets aroun d low mass stars dRavmond et al.1120071 : 
lLissaueil[20 07). Indeed, the cases where we followed tra¬ 
ditional assumptions by not increasing the surface den¬ 
sity of solids at the water line, we found dry planets in 
the habitable zones for stars with masses <0.6 Mq as in 
previous work. However, in those cases where the solid 
surface density was increased at the water line (the IC 
cases), we find that dry habitable zone planets may be 
limited to only very low mass disks or stars ^0.2 Mq or 
less. 

The results of these types of models depend sensitively 
on the choice of initial conditions used in these simula¬ 
tions, such as the water line location or surface density 
of planetesimals. In order to avoid exploring a large pa¬ 


rameter space, we limited ourselves to only one surface 
density profile for the building blocks of the planets and 
linearly scaled the surface density of the planetary build¬ 
ing blocks with the central mass of the star. This linear 
scaling is consistent with the pr otop l anetary di sk mass 
scaling relationship reported by lAndrews et al.l (1 2013 1, 
though significant variations among disks demonstrate 
such scaling should be used as a guide a nd not a defini te 
rule. This issue was also discussed by iRavmond et al.1 
(2007;. and like those authors, we suggest that these 
results be taken as a statistical guide for inferring the 
properties of planets around stars of different masses. 

While we did not model a range of disk masses around a 
single stellar mass, we can use the results presented here 
to draw some inferences about how planetary properties 
may vary for disks with masses different from those con¬ 
sidered here. One trend to note is that the ranges in the 
properties of the habitable zone planets also scale with 
mass of the central star. Figure 5 shows how the masses 
and water mass fractions of the planets in the habitable 
zones of the respective stars vary in the runs presented 
here. For the lowest mass stars considered, 0.2 Mq and 
0.4Mq, we find that the masses of the planets in the hab¬ 
itable zone are tightly clustered at values that are just 
slightly in excess of the mass of the embryos used in those 
simulations, with variations of only ~10-20%. That is, 
these planets generally represent single embryos which 
accreted a number of planetesimals and avoided impacts 
with other embryos. This was true in both the CC cases 
and IC cases, as the planets that formed in the habit¬ 
able zones around these stars were sufficiently far from 
the location of the water line that the increase in surface 
density in the further reaches of the disk did not have a 
dramatic effect on their growth. The more massive stars, 
however, have habitable zone planets whose masses range 
by factors of 2-20. Likewise, the water mass fraction of 
habitable zone planets cluster at minimal values for the 
low mass stars (all have WMF=10 -5 in the 0.2 Mq case), 
but the range increases with increasing mass. The range 
in WMF is actually greater in the CC cases than the 
IC cases for stars >O.6M 0 , likely because the fraction 
of water-bearing building blocks is significantly greater 
in the IC cases, meaning individual accretion events can 
change the water inventory of a planet more significantly 
than in the CC cases. 

The increase in the ranges in planetary mass with in¬ 
creasing stellar mass seen in the simulations is a product 
of the chaotic nature of planetary accretion caused by 
the gravitational encounters of a number of randomly 
placed embryos. In the low stellar mass cases, gravita¬ 
tional encounters among planetary embryos are less fre¬ 
quent and less intense due to the fact that mass is less 
available and less concentrated in these low mass disks. 
Thus planetesimal-embryo or embryo-embryo collisions 
and scattering will not occur as frequently as in disks 
with higher masses, allowing for more orderly growth. 
Thus, the range in properties of the planets seen in these 
models is likely a product of the disk properties and less 
so the properties of the star. This suggests that planets 
that form from more massive disks are likely to display a 
wider range of characteristics than those that form from 
lower mass disks. Further, if this extrapolation is cor¬ 
rect, terrestrial planets which form from low mass disks 
are likely to serve as fossil records of the the properties 
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of the protoplanetary disk-the low levels of gravitational 
scattering among planetary building blocks mean that 
planets which form in such disks would grow from plan- 
etesimals and embryos that formed locally. 

An important caveat to this work is that the ef¬ 
fects of more massive planets, like gas or ice giants 
and super-Earths, have been neglected. In our own So¬ 
lar System, different orbital properties of Jupiter and 
Saturn would have led to significantly different accre¬ 
tion histories for the inner planets, differences which 
go beyond those at tributed to the stochastic nature of 
planetary accretion (lO’Brien et a l. 2006; Raymond et al.1 

120091 [Fischer fc Cieslall20l4lQuintana fe Lissaueill2014f) . 

In addition, migration of gas giants could play a crit¬ 
ical role in determining how ter restrial planets accrete 
in a young plan etary system (|Ravmond et, all 120061 : 
iWalsh et al.ll2f)l ill . Thus such planets, if present, could 
lead to different accretion outcomes and different rela¬ 
tionships between the habitable zone planets and the 
stellar mass than found here. Again, we neglected 
such effects in order to focus on the consequences of 
the increased surface density at the water line. How¬ 
ever, gas giants are not present in all systems; in fact, 
their occurrence rate may be as low as 10%, meaning 
most systems a re likely to be born without any such 
planets present feasper et al.ll2007t IHoward et al.ll2012t 


Mortier^_aL 1120121 iZechmeister etal .1120131 iBiller et al.1 


201 . 'll IDressing fc CharbonneaiJl201,Jfl . 

Super-Earths and ice giants are present around stars 
at higher frequencies than gaseous gi ants, t hough still 
are not present in all systems (IHoward et all [20101 

IDressing fc Charbonneaull2013llFressin et al .1120131) . The 

formation of these objects remains uncertain. It is possi¬ 
ble that such objects form through accretion of solids, as 
outlined here, coupled with radial migration through the 
gaseo us disk ([Raymond et alJl2008llRavmond fc Cossoul 
12014) . Indeed, depending on the rate of migration, the 
inward movement of accreting bodies could play a impor¬ 
tant role in setting the structure of the planetary system 
and chemical compo sition of the planets that form (e.g. 
IQgihara fc Idal 120 091. Indeed, the effects of planetary 
migration on planetary accretion is likely variable due to 
uncertainties on drift rates and will likely vary depend¬ 
ing on turbulence levels and the mass of the disk at the 
time accretion begins (lChambersll2009l) . 

Others have suggested that super-Earths formed in situ 
through gravitational accumulation of planetesimals an d 
embryo s as described here (IChiang fc Laughlinl l2fll.1l '). 
Indeed, Kennedy et al.1 1)200611 explored the formation of 
such planets using analytic models for accretion in the 
context of the framework described here; that is, such 
planets were thought to possibly form outside the water 
line of a protoplanetary disk where the surface density of 
solids was increased by a factor of 2 or 3 due to the con¬ 
densation of water ice as considered here. We have shown 
that one would expect more massive planets to form 
when such a surface density increase is accounted for, 
though to get very massive planets would likely require 
much more massive disks than assumed here. Again, this 
should be the focus of future work. 

Thus the detailed results here are most applicable to 
those planetary systems which form from relatively low 
mass disks and where massive planets are absent. De¬ 
spite the additional effects that could be considered in 


the future, this work has shown that the increase in 
surface density and high water-content of bodies ex¬ 
pected beyond the water line of a protoplanetary disk 
can have an important effect on the accretion of ter¬ 
restrial planets. Specifically, the increase in mass be¬ 
yond the water line will lead to more gravitational scat¬ 
tering than in those cases where bodies are treated as 
carbonaceous chondrite analogs, leading to larger and 
more volatile-rich planets than seen in previous stud¬ 
ies. The water-rich planetary building blocks consid¬ 
ered here are likely the norm in planetary systems 
which contained little to no short-lived radionuclides. 
Even in our own solar system, planetary building blocks 
which originated beyond the water line were likely more 
water-rich than the carbonaceous chondrites we see to¬ 
day as aqueous alteration, fluid flow, and water loss 
would have taken millions to tens of millions of years to 


occur ( Grimm & McSweed 1989f 1 Cohen & Coker 

2000 ; 

Young et al. 20031 Castillo-Rogez & McCord 

20 m 


Thus during the period of planetary assembly, these bod¬ 
ies likely looked different than previously assumed, allow¬ 
ing for water-rich bodies to be more common than had 
been realized. The effects of having more water-rich par¬ 
ent bodies for the carbonaceous chondrites in our own 
Solar System need to be investigated. 

In considering the sui te of results presented h ere, along 
with those described in lRavmond et all (120071) . the gen¬ 
eral trends that emerge are interesting to consider as we 
collect more data on the properties of extrasolar plan¬ 
ets and their planetary systems. As discussed above, 
these simulations predict two major trends to develop: 
planetary systems will contain greater numbers of lower 
mass planets in small orbital spacing around low mass 
stars while higher mass stars will be home to fewer, more 
massive planets. In addition, within the simulations for 
a single stellar mass, the predicted planetary properties 
follow general trends showing how water mass fraction 
would vary with semi-major axis, as shown in Figure 6. 
Thus, even though every planetary system that forms in 
these studies is different, due to the different assumed po¬ 
sitions of the initial building blocks, the resulting suite 
of planets define a radial trend. Such predictions pro¬ 
vide an opportunity to test N-body accretion studies, 
or the choices of initial conditions used by such stud¬ 
ies, against actual data. That is, as shown in previ¬ 
ous work, the outcome of individual N-body simulations 
are highly dependent on the initial conditions, meaning 
that individual runs are unlikely to reproduce all of the 
properties of any given planetary system. However, tak¬ 
ing the collection of simulations together and comparing 
them to the suite of exoplanetary systems, specifically 
the IFO/roek ratio as determined by mass-radius rela¬ 
tionships le.g.lSeager et al.ll2007l [Rogers fc Seageill2010l 
iSwift et al.ll2fll2D . could allow us to test our ideas of this 
stage of planet formation in greater detail. Should the 
trends predicted here be consistent with the planetary 
systems that are detected, then this would give support 
that our choices for the initial conditions in these N- 
body simulations are accurately describing those from 
which planetary systems develop. Should we see trends 
that are inconsistent with these predictions, that would 
suggest that the way in which we initialize our N-body 
simulations may not accurately describe the conditions 
under which planetary accretion begins, requiring other 
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possibilities to be considered in detail. In fact, trends 
such as these may allow us to infer the location of the 
water line in protoplanetary disks as shown in Figure 6— 
the water line was located at 0.95 AU in these simulations 
and one readily sees that this point marks a difference 
between planets that have nearly uniform water mass 
fractions (outside the water line, and thus the plateau 
in this plot) and those that have decreasing water con¬ 
tents as one gets closer to the star (inside the water line). 
Alternatively, disagreement in data and the model pre¬ 
dictions could indicate that the processes ignored in N- 
body simulations, such as inelastic collisions and perfect 
mergers or disk-driven migration, are actually important 
and must be considered in accretion studies. Thus the 
growing collection of constraints on the numbers, masses, 


and orbital properties of extrasolar planets now allows us 
to explicitly test our models for planet formation. Such 
efforts will help us to better understand the origin of our 
own Solar System. 
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Table 1: Simulation Parameters 


M* 

(Mq) 

Initial 

Range 

(AU) 

Water 

Line 

(AU) 

M e mb 

(M 0 ) 

M p lan 

(M 0 ) 

N cc 

emb 

N CC 

plan 

N ic h 

emb 

Ad, c 

plan 

Habitable 

Zone 

(AU) 

1.0 

0.5-4.0 

1.3 

0.05 

2.5xl0" 3 

50 

974 

88 

1875 

0.8-1.5 

0.8 

0.2-3.2 

1.1 

0.04 

2xl0 -3 

43 

827 

72 

1528 

0.39-0.74 

0.6 

0.1-2.4 

0.9 

0.03 

1.5xl0 -3 

33 

622 

54 

1122 

0.20-0.37 

0.4 

0.05-1.6 

0.8 

0.02 

10“ 3 

23 

420 

34 

672 

0.10-0.19 

0.2 

0.03-0.8 

0.5 

0.01 

5xl0 -4 

11 

196 

16 

291 

0.05-0.10 


Embryo mass; M p i an : Planetesimal mass; AT em ^: Number of embryos, N p i an : Number of planetesimals. Habitable Zone estimates 
taken from I Raymond et alJ (|2007l ) 



Time [yrs] 


Fig. 1. — Location of the water line, defined by T= 140 K, as a function of time for 0.2, 0.4, 0.6, 0.8, and 1.0 Mq stars based on the 
models of ISiess et al.l (|2000n . 
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Fig. 2.— Series of snapshots showing the accretion history of one of our OAMq CC cases. Symbol size scales with the mass of the bodies 
shown, while colors indicate the corresponding water mass fraction. 
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log(Water Mass Fraction) 
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Fig. 3.— Series of snapshots showing the accretion history of one of our 0.4 Mq IC cases. Symbols and colors are as shown in Figure 2. 
An important difference to note here is that the water mass fraction extends to higher values here than in the CC cases. Also note that the 
surface density of bodies beyond the water line in this simulation is greater than the CC case. This leads to more dynamical excitement of 
bodies in the outer disk and, thus, more radial mixing of planetary building blocks. 
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Fig. 4. — Final architectures of all simulations done in this study, in order of increasing stellar mass. Vertical lines indicate the location 
of the Habitable Zone as given in Table 1. As in Figures 2 and 3, symbol size scales with planetary mass to the ^ power, while symbol 
color indicates the water mass fraction of the planet. Traditional, CC cases are shown in the bottom of each simulation, while the more 
water-rich, IC cases are shown above them. The Solar System planets are shown for reference in each panel. 
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Star Mass [Solar Mass] 


Fig. 5. — Properties of only those planets which reside in the Habitable Zone for each star at the end of the simulations as a function of 
stellar mass. Top panel shows the masses of the Habitable Zone planets, and how such planets exhibit greater ranges in values for higher 
mass stars. This is due to the more extensive gravitational interactions that occur in the higher mass disks which accompany such stars, 
and thus the more chaotic nature of this evolution. Growth around lower mass stars is much more orderly. The lower panel shows the 
water mass fractions for the Habitable Zone planets in this simulation, showing that the low mass stars, due to the low masses of the disks 
around them, do not see significant inward scattering of water-rich bodies during planetary accretion, while inward scattering is occurs 
more readily in the more massive disks (either the IC cases or higher mass stars). 
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Fig. 6.— Cumulative distribution of planetary properties for the 0 .6M@ runs performed in this study, showing how planetary properties 
vary as a function of semi-major axis. While every simulation produced different planetary systems, a trend that develops in these 
simulations is the radial gradient in the water mass fraction of the planets that arise from radial diffusion of water-bearing planetesimals 
due to gravitational interactions beyond the water line. Looking for trends like these in the collection of exoplanetary systems would serve 
as important tests for our models and ideas on the final stages of planetary accretion. 



